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Chapter  1 


Introduction 


The  IC  market  is  characterized  by  an  ever  increasing  level  of  integration  com¬ 
plexity.  Today  complete  systems  that  previously  occupied  one  or  more  boards,  are 
integrated  on  a  few  chips  or  even  on  one  single  multi-million  transistor  chip  -  a  so 
called  System-  on-  Chip  (SoC). 

Research  interest  is  moving  in  the  direction  of  system  synthesis  where  an  object- 
oriented  system  level  specification  is  translated  into  a  hardware-software  co-architecture 
with  high  level  specification  for  both  the  hardware,  the  software  and  the  interfaces. 
In  addition,  reuse  and  platform  based  design  methodologies  are  being  developed  to 
further  reduce  the  design  effort  for  complex  systems. 

Millimeter-wave  circuits  are  becoming  ever  more  commercially  and  militarily  vi¬ 
able  and,  coupled  with  large  scale  production,  design  practices  must  evolve  to  be 
more  sophisticated  in  order  to  handle  more  complex  relationship  between  constituent 
elements.  Also,  as  the  frequency  of  radio  frequency  (RF)  circuits  extends  beyond  a 
gigahertz  to  tens  and  hundreds  of  gigahertz,  wavelengths  become  large  with  respect 
to  device  and  circuit  dimension  and  the  three-dimensional  EM  environment  becomes 
more  signihcant. 

Microwave  and  millimeter-wave  systems  have  generally  been  developed  from  a 
circuit  perspective  with  the  effect  of  the  electromagnetic  (EM)  environment  modelled 
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using  lumped  elements  or  iV-port  scattering  parameters.  The  recent  development  of 
the  local  reference  terminal  concept  coupled  with  steady-state  and  transient  analysis 
nsing  state  variables  allows  the  incorporation  of  nnrestrained  EM  modelling  of  mi¬ 
crowave  strnctnres  in  a  circnit  simnlator.  A  strategy  implementing  global  modelling 
of  electrically  large  microwave  systems  using  the  circnit  abstraction  is  presented. 


1.1  Motivation  for  and  Objective  of  this  Study 

One  fnndamental  issne  that  arises  in  modelling  electrically  large  systems  in  a  cir¬ 
cnit  representation  is  the  assignment  of  system  gronnd.  Cnrrently,  circnit  simnlators 
use  a  nodal  approach  in  which  voltages  are  assigned  to  terminals  and  each  of  these 
voltages  is  referred  to  a  common  reference  point,  which  is  commonly  called  ‘ground’. 
In  a  spatially  distribnted  system,  a  common  reference  point  cannot  always  be  dehned 
as  the  (electrically  signihcant)  spatial  separation  of  a  node  and  its  reference  point 
cannot  be  tolerated.  If  the  separation  is  an  appreciable  fraction  of  a  wavelength,  it  is 
not  possible  to  nniqnely  dehne  voltage. 

Transmission  line  effects  frequently  lead  to  problems  like  overshoot,  nndershoot 
and  ringing.  Board  interconnects  limit  performance  by  introdncing  delay  into  the 
timing  path.  Crosstalk  and  EMI  are  beginning  to  present  formidable  problems  in 
what  traditionally  have  been  considered  acceptable  board  designs. 

EM  models,  relying  on  the  global  reference  terminal  concept  cannot  accurately 
treat  large  multi-port  structures.  This  thesis  presents  a  tool  flow  for  model  develop¬ 
ment  based  on  a  new  concept  of  local  reference  nodes  (LRN’s).  The  LRN  concept 
enables  the  conversion  of  a  port  -  based  model  to  a  terminal-  based  model,  which 
can  be  used  in  a  general  pnrpose  circnit  simulator  after  modihcations  to  the  way  the 
nodal  admittance  matrix  is  treated. 

There  is  a  lot  of  research  being  carried  ont  for  the  simnlation  of  qnasi-optical 
grid  ampliher  strnctnres.  A  Method  of  Moment  simnlator  was  created  by  Nnteson 
[9]  which  generates  a  mnlti-port  impedance  model  of  the  passive  section  of  a  grid 
structure.  Summers  [8]  describes  a  tool  flow  for  linking  the  passive  strnctnres  and 
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the  active  devices  in  a  non-linear  circuit  solver.  Since  the  system  consists  of  spatially 
distributed  components,  we  need  a  system  which  can  identify  the  local  reference 
terminals. 

This  thesis  will  outline  the  creation  of  a  tool  set  with  the  capability  to  simu¬ 
late  spatial  power  combining  systems.  This  thesis  documents  the  development  of  a 
tool  flow  for  integrated  steady  state  nonlinear  circuit  analysis  and  electromagnetic 
analysis  of  a  quasi-optical  grid  ampliher  system.  The  analysis  incorporates  surface 
modes,  nonuniform  excitation,  and  full  nonlinear  effects.  The  work  is  verihed  using 
measurement  of  a  5x5  grid  ampliher. 


1.2  Thesis  Overview 

This  thesis  describes  the  integration  of  a  circuit  and  held  simulator,  which  will 
ultimately  be  used  to  simulate  a  grid  ampliher. 

Chapter  2  is  a  review  of  literature  relevant  to  the  grid  ampliher  systems.  The 
review  includes  modeling  of  quasi-optical  arrays,  local  reference  terminal  concept  for 
simulating  spatially  distributed  circuits  and  the  need  for  a  technology  hie  to  describe 
the  process  in  a  multilayered  system. 

Chapter  3  is  an  overview  of  the  Java  philosophy  and  explains  the  programming 
interface. 

Chapter  3  is  an  overview  of  the  program  how  required  to  complete  the  simula¬ 
tion.  The  overview  includes  a  discussion  about  each  program  used  and  a  detailed 
description  of  the  interfaces  between  the  programs  are  presented. 

Chapter  4  is  a  case  study  on  modeling  a  grid  ampliher.  A  detailed  description 
of  the  model  development  process  as  well  as  results  from  the  completed  model  are 
presented 

Chapter  5  presents  conclusions  derived  from  this  project  and  areas  that  need 
further  exploration  in  the  future. 
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Chapter  2 


Literature  Review 


2.1  Background 

Quasi-optical  power  combining  techniques  are  a  means  of  combing  the  power  out¬ 
put  of  numerous  solid  state  devices  generally  in  a  free  space  environment.  In  [6] 
a  review  of  the  various  techniques  and  systems  is  presented.  Mink  [3]  conducted 
a  theoretical  investigation  of  an  array  of  hlamentary  current  sources  radiating  into 
a  plano-concave  open  resonator.  Mink  [3]  discuses  the  possibility  of  obtaining  high 
combining  efficiency  for  systems  several  wavelengths  in  dimension.  The  held  was  fur¬ 
ther  advanced  by  the  construction  of  a  grid  oscillator  [4]  and  grid  ampliher  [5].  A 
diagram  of  the  general  system  is  shown  in  Fig. 2.1  The  systems  represented  a  proof  of 
concept. 

The  number  of  spatial  power  combining  systems  has  increased  greatly  over  the 
years,  but  the  design  of  the  systems  has  largely  been  empirical  or  through  the  use  of 
simplihed  models.  The  task  of  characterizing  the  structures  has  proved  to  be  difficult 
because  the  systems  are  usually  several  wavelengths  in  dimension.  Some  signihcant 
work  has  been  done  in  the  area  of  planar  grid  structures.  Early  models  assumed  the 
grid  was  inhnitely  large,  but  recent  developments  [9],  have  removed  this  assumption 
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OUTPUT  HORN 


Figure  2.1:  Planar  Grid  System. 

from  the  model.  The  literature  review  discusses  modeling  techniques  for  quasi-optical 
grid  amplihers  and  oscillators. 

The  power  combining  system  is  a  way  of  efficiently  coupling  the  energy  via  free 
space  instead  of  using  transmission  lines.  If  you  use  transmission  lines,  there  is  an 
inherent  upper  limit  on  the  number  of  elements  that  can  be  combined.  It  hnds 
application  in  X-band  (8  Ghz)  and  Ka  -  band  (34Ghz)  combiners. 


2.2  Modeling  of  Quasi-Optical  System 

The  grid  oscillator  presented  in  [4]  was  one  of  the  early  demonstrations  of  quasi- 
optical  power  combining.  The  system  was  designed  using  a  unit  cell  model  that  is 
presented  in  [9].  In  [5]  Kim  presents  a  grid  ampliher.The  equivalent  circuit  and  unit 
cell  layout  are  shown  in  Fig.2.2  and  Fig.  2.3. 

In  [9]  the  unit  cell  model  was  extended  to  improve  its  performance.  The  surface 
current  distribution  was  approximated  using  the  method  of  moments  technique  and 
the  EMF  method  was  used  to  calculate  the  impedance  of  the  grid  leads.  The  sub¬ 
strates  and  free  space  regions  in  the  systems  are  represented  by  transmission  lines. 
The  polarizers  are  modeled  by  reactive  elements. 

Research  at  North  Garolina  State  University  led  to  the  development  of  a  Method 
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Figure  2.2:  Unit  cell  with  gate  bias  lines. 


Figure  2.3:  Equivalent  circuit  model. 
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of  Moments  simulator  [9]  which  no  longer  relied  on  the  unit  cell  approximation  in 
characterizing  planar  grid  structures.  The  signihcance  of  this  work  is  that  surface 
waves,  nonuniform  excitation,  and  gain  and  phase  variations  of  the  unit  amplihers 
can  be  considered.  Finite  grids  can  now  be  rigorously  specified,  with  the  effects  of 
mutual  coupling  between  cells  and  variations  in  magnitude  and  phase  of  the  feed 
sources  included.  The  simulator  also  included  the  effects  of  quasi-optical  components 
far  away  from  the  grid,  such  as  lenses.  The  simulator  computed  an  impedance  ma¬ 
trix  for  the  grid  structure,  which  can  be  used  in  a  microwave  circuit  simulator.  An 
implementation  is  presented  in  [10].  The  gain  of  an  ampliher  grid  and  an  analysis  of 
frequency  and  phase  synchronization  in  oscillator  arrays  are  presented. 


2.3  Local  Reference  Terminal 

Circuit  theory  has  evolved  to  include  a  common  reference  terminal  to  which  all 
voltages  in  a  circuit  are  referred.  However,  it  is  generally  not  feasible  to  dehne  voltages 
or  a  single  reference  point  in  a  spatially  distributed  system.  In  [1]  ,  Steer  describes 
a  framework  to  simulate  distributed  strucutres.  Microstrip  networks  are  examples  of 
distributed  structures  for  which  reasonable  approximations  have  been  made  so  that 
they  can  be  treated  as  elements  in  conventional  circuit  simulators.  A  micro-strip 
transmission-line  segment  is  shown  in  Fig.  2.2  with  the  conventionally  accepted  deh- 
nition  shown  for  the  node  voltages  at  the  ends  of  the  line.  In  the  common  approach, 
the  voltages  Vi  and  V2  are  determined  as  the  integral  of  the  electric  field  from  the 
strip  to  the  ground  plane  using  the  shortest  path  -  defined  with  vertical  arrows.  How¬ 
ever,  an  absurd  ‘value’  for  V2  would  be  obtained  if  the  integral  were  to  be  performed 
along  the  dashed  line.  Thus,  the  common  approach  treats  the  ground  plane  itself 
as  a  global  reference  point  even  though  this  implies  and  requires  that  charges  can 
instantaneously  redistributed  on  the  ground  plane.  While  this  has  proven  to  be  a 
successful  modelling  approach  for  microwave  microstrip  circuits  that  are  generally 
arranged  as  cascaded  stages,  limitations  are  becoming  evident  in  the  modelling  of 
strip-like  distribution  networks  in  high  speed  digital  systems. 
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Figure  2.4:  A  two  port  micro-strip  line. 

2.4  Technology  file 

In  multi-layered  structures,  we  need  to  define  the  layers  that  are  used  in  the  sys¬ 
tem.  The  technology  file  contains  process  specific  parameters  such  as  layer  thicknesses 
and  the  sheet  resistance  of  the  various  layers.  In  addition  to  this  it  also  provides  the 
mesher  tool  with  information  regarding  the  ordering  of  layers.  A  generic  CMOS  pro¬ 
cess  is  shown  in  Figure  2.5.  In  Figure  2.6  a  sample  technology  file  for  the  given  process 
is  shown. 

In  conventional  layout  tools,  the  technology  file  contains  information  pertaining 
to  lambda-based  design  rules.  They  do  not  store  the  material  properties  (dielectric 
constant),  (loss  tangent),  (sheet  resistivity),  (conductivity),  which  are  essential  in 
characterizing  multi-layered  structures.  The  main  advantage  to  having  this  informa¬ 
tion  external  is  that,  we  can  change  our  process  parameters  without  having  to  redo 
our  layout. 

In  addition  to  this,  the  information  pertaining  to  the  order  of  layers  is  also  embed¬ 
ded  in  the  file.  This  helps  us  render  a  two-dimensional  layout  to  a  three-dimensional 
layout. 

The  output  of  the  program  is  a  text  file,  which  is  human  readable.  The  informa¬ 
tion  available  in  the  file  can  be  parsed  to  generate  inputs  for  other  layout  tools  like 
Cadence,  LinkCad,  L-Edit  etc. 
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Figure  2.5:  A  sample  process  parameters. 
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//TECHNOLOGY  FILE 

II  METAL  LAYERS  :  2 
II DIELEOTRIO  LAYERS:  1 

technology  cmos 
UNITS  microns 
metal 

//name  zmln  thickness  reslstlvlty(ohms,cm)  desc 

layer  0  0,00  0,35  0,2  //G  substrate 

layer!  0,7  0,2  l,56E-4  II M 

endmetal 

dielectric 

//name  zmln  thickness  dielectric  constant 

layer  2  0,35  0,35  3,9 

enddlelectric 
endtechnology 


Figure  2.6:  Sample  technology  file. 
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Chapter  3 


Java  Philosophy 


3.1  Introduction 

Java  is  a  system-independent  object  oriented  programming  language.  The  pro¬ 
gram  architecture  is  based  on  the  MVC  (Model-View  -Controller  architecture).  The 
goal  of  the  MVC  design  pattern  shown  in  Figure  3.1  is  to  separate  the  application 
object  (model)  from  the  way  it  is  represented  to  the  user  (view)  from  the  way  in 
which  the  user  controls  it  (controller). 

The  model  object  knows  about  all  the  data  that  need  to  be  displayed.  It  also  knows 
about  all  the  operations  that  can  be  applied  to  transform  that  object.  However,  it 
knows  nothing  whatever  about  the  GUI,  the  manner  in  which  the  data  are  to  be 
displayed,  nor  the  GUI  actions  that  are  used  to  manipulate  the  data.  The  data 
are  accessed  and  manipulated  through  methods  that  are  independent  of  the  GUI.  A 
simple  example  of  a  model  would  be  a  clock  object.  It  has  intrinsic  behaviour  whereby 
it  keeps  track  of  time  by  updating  an  internal  record  of  the  time  ever  second.  The 
object  would  provide  methods  which  allow  view  objects  to  query  the  current  time.  It 
would  also  provide  methods  to  allow  a  controller  object  to  set  the  current  time. 

The  view  object  refers  to  the  model.  It  uses  the  query  methods  of  the  model  to 
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obtain  data  from  the  model  and  then  displays  the  information.  The  display  can  take 
any  form,  in  the  clock  example,  one  view  object  could  display  the  time  as  an  analogue 
clock,  another  could  show  it  as  a  digital  clock.  The  different  displays  would  have  no 
bearing  whatsoever  on  the  intrinsic  behavior  of  the  clock. 

The  controller  object  knows  about  the  physical  means  by  which  users  manipulate 
data  within  the  model.  In  a  GUI  for  example,  the  controller  object  would  receive 
mouse  clicks  or  keyboard  input  which  it  would  translate  into  the  manipulator  method 
which  the  model  understands.  For  example,  the  clock  could  be  reset  directly  by  typing 
the  current  time  into  the  digital  clock  display.  The  controller  object  associated  with 
the  view  would  know  that  a  new  time  had  been  entered  and  would  call  the  relevant 
‘SetTime’  method  for  the  model  object. 


3.2  Event  Delegation  Model 

The  delegation  event  model  supports  a  clean  separation  between  the  core  program 
and  its  use  interface.  The  delegation  pattern  allows  robust  event  handling  that  is  less 
error  prone  due  to  strong  compile-  time  checking. 

The  delegation  model  consists  of  three  separate  objects  to  deal  with  event  han¬ 
dling:  event  sources,  events,  and  listeners.  Event  sources  are  any  of  the  user  interface 
components  such  as  buttons,  text  helds.  Event  sources  generate  and  hre  events,  and 
the  events  propagate  until  they  are  hnally  attended  to  by  their  listeners. 


3.3  Class  Structure 

Model  Class  The  model  class  is  GIF  hie  parser  and  the  Technology  hie  generator. 
View  Class  The  view  class  is  composed  of  the  diherent  panels  (ex.  metalPanel, 
dielectricPanel,  viaPanel,  maiuMenu) 

Controller  Class  The  controller  class  is  embedded  in  the  view  class.  Each  indi¬ 
vidual  component  has  its  own  listener  and  adapter  class. 
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3.4  Summary 

We  have  chosen  Java  as  the  programming  language  as  it  simplifies  the  program¬ 
ming  effort  by  providing  a  single  file  which  can  run  on  different  platforms.  Since  Java 
is  based  on  the  event-delegation  model,  this  allows  us  to  separately  handle  the  events 
generated  by  individual  objects  in  the  program. 
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Chapter  4 


Tool  Flow 


4.1  Introduction 

Figure  4.1  is  a  typical  spatial  power  combining  system.  The  system  contains  horns 
for  input  and  output  and  components  which  are  spread  over  several  wavelengths.  The 
active  surface  contains  numerous  nonlinear  devices  connected  to  input  and  output 
radiators. 

The  system  presents  some  very  challenging  problems:  components  are  spatially 
distributed,  power  is  radiated  through  the  system,  and  nonlinear  devices  are  present 
in  the  system.  The  hrst  two  problems  dictate  that  electromagnetic  models  be  used, 
however  the  nonlinear  devices  are  naturally  handled  in  a  circuit  simulator.  The 
solution  is  to  divide  the  two  problems  up  as  shown  in  Figure  4.2  . 

The  spatially  distributed  elements  are  handled  through  electromagnetic  simula¬ 
tion,  and  the  nonlinear  devices  are  handled  by  common  circuit  simulation  techniques. 
The  interfaces  between  the  two  worlds  must  be  clearly  dehned.  The  circuit  simula¬ 
tor  requires  the  impedance  network  seen  by  the  active  devices.  The  simulator  also 
requires  any  excitation  currents  or  voltages  generated  by  input  helds.  The  excitation 
should  be  normalized  to  limit  the  amount  of  electromagnetic  simulations  required. 
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The  electromagnetic  simulator  will  require  the  voltages  at  the  ports  of  the  active 
devices  to  calculate  the  current  distribution  and  radiated  helds. 


nillPIITPOIARI/FR 


Figure  4.1:  Typical  Grid  Ampliher  System. 


Figure  4.2:  Division  of  system  in  fields  and  circuits. 
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4.2  Local  Reference 

Once  the  division  of  the  system  is  done,  we  need  to  identify  the  local  reference 
terminals  and  the  local  reference  groups.  The  port-based  admittance  matrix  is  nodal 
in  nature  as  illustrated  in  Figure  4.3.  The  local  reference  terminals  are  marked  as 
there  is  no  global  ground.  The  implementation  of  the  local  reference  terminal  is 
discussed  in  more  detail  in  [2]. 


LOCAL  REFERENCE 


Figure  4.3:  Unit  cell  with  port  voltages  and  local  reference  terminals  identihed. 
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In  our  simulation  we  assume  each  unit  cell  with  the  active  device  to  form  the  local 
reference  group  as  shown  in  Figure  4.4. 


Figure  4.4:  Local  Reference  Groups  shown  on  a  4x4  grid  ampliher. 
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4.3  Tool  Flow 


Figure  4.5  shows  the  complete  block  diagram  of  the  simulation  environment  which 
is  capable  of  handling  the  distributed  effects  in  high  speed  circuits.  Off  the  shelf  tools 
were  taken  and  modihed  to  capture  the  distributed  effects.  The  main  contribution  to 
the  tool  flow  was  adding  support  to  capture  the  geometrical  and  electrical  properties 
of  a  layout  and  handling  the  interface  between  the  programs. 


Figure  4.5:  Block  Diagram  of  Tool  Flow. 


4.3.1  Simulating  High  Speed  Circuits 

Figure  4.6  shows  the  steps  involved  in  preparing  a  generic  layout  for  EM  sim¬ 
ulation.  The  system  under  observation  is  drawn  using  Electric.  The  geometrical 
information  is  automatically  taped  out  as  a  GDS-II  hie.  The  additional  information 
included  in  the  GDS-II  hie  is  the  local  reference  terminals  that  are  manually  identihed 
by  the  user. 

The  primary  input  to  EMPDK  is  the  circuit  layout  description  in  GIF  or  GDS-II 
format.  This  hie  is  parsed  and  the  user  is  presented  with  the  layers  in  the  drawing. 
The  user  can  then  manually  select  the  layers  of  interest  and  provide  the  electrical 
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parameters  required  for  EM  -  simulation.  This  stage  prepares  the  layout  hie  for  EM 
-  simulation. 

The  output  hies  generated  by  Electric  and  EMPDK  completely  describe  a  dis¬ 
tributed  structure.  This  information  is  then  fed  to  the  held  solver  which  provides 
us  with  the  port-based  parameters.  The  port-based  parameters  characterize  the  dis¬ 
tribute  structure  and  this  information  is  used  to  simulate  circuits  with  active  devices. 


Add  process  information 
for  each  layer 


Generate 

technology 

file 
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4.3.2  Schematic  Capture  Flow 

The  tool  flow  presented  in  Figure  4.5  can  also  be  used  to  study  simple  circuits 
with  microwave  elements  directly  from  their  schematic.  fREEDA  has  a  good  set  of 
microwave  elements.  But  it  is  a  circuit  netlist  interpreter.  To  make  the  simulation 
environment  more  user  friendly,  we  need  a  graphical  user  interface.  This  allows  the 
user  to  draw  the  schematic  and  automatically  generate  the  netlist  for  simulation.  At 
this  stage,  the  user  has  to  manually  enter  the  output  statements  and  the  analysis 
type  required  in  the  simulation.  The  flow  is  explained  in  Figure  4.7. 
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4.4  EMPDK 

4.4.1  Introduction 


Figure  4.8:  EMPDK  vl.O 

EM-Physical  Design  Kit  is  a  Java  program  which  takes  generic  GDS-II  or  GIF  hie 
as  its  input  and  produces  an  output  GIF  hie  and  an  associated  technology  hie.  This 
progam  was  developed  as  part  of  the  work  for  this  thesis.  It  provides  an  interface 
between  the  layout  tool  and  full-wave  electromagetics  simulator. 

The  technology  hie  gives  us  information  regarding 


•  the  material  properties  of  the  layers  in  the  process 
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•  the  ordering  of  the  layers,  which  helps  us  render  the  layout  to  a  three-dimensional 
view. 

•  provides  thickness  information  of  each  layer. 

The  geometrical  information  of  the  on-chip  placement  is  obtained  from  the  GIF  hie. 
The  technology  hie  gives  the  mapping  between  the  diherent  layers  in  the  GIF  hie  and 
the  generic  layers  that  any  layout  tool  can  interpret.  The  output  of  EMPDK  is  a 
human  readable  text  hie. 


4.4.2  Creating  a  technology  file 

//TECHNOLOGY  FILE 

II  METAL  LAYERS  :  2 
II DIELEOTRIG  LAYERS:  1 

technology  cmos 
UNITS  microns 
metal 

//name  zmln  thickness  reslstlvlty(ohms,cm)  desc 

layer  0  0,00  0,35  0,2  //G  substrate 

layer  1  0,7  0,2  l,56E-4  II M 

endmetal 

dielectric 

//name  zmln  thickness  dielectric  constant 

layer  2  0,35  0,35  3,9 

enddlelectric 
endtechnology 


Figure  4.9:  Sample  technology  hie 

A  sample  technology  hie  is  shown  is  Fig  4. 9. The  ‘layer’  statements  dehne  the 
various  dielectric  and  substrate  layers.  A  conducting  substrate  layer  is  specihed  by 
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stating  its  resistivity  in  ohm-cm.  A  dielectric  layer  can  be  modeled  as  an  equiva¬ 
lent  resistivity  according  to  the  dielectric  loss  tangent.  The  thickness  is  specihed  in 
microns.  The  relative  permitivity  is  unitless.  Notice  that  the  substrate  layers  are 
dehned  from  bottom  up.  In  other  words,  the  order  is  important.  Thus  the  bulk  layer 
is  listed  hrst,  followed  by  the  epi,  and  hnally  followed  by  the  oxide  layer. 

The  conducting  metal  layers  are  dehned  in  the  ‘metal’  sections.  Each  metal  layer 
belongs  to  a  specihc  layer  specihed  with  the  layer  identihcation.  Most  metal  layers 
will  be  dehned  in  non-conducting  oxide  layers.  The  sheet  resistance  of  the  metal  layers 
is  specihed  in  milli-Ohms/Square  along  with  the  metal  thickness.  The  z-coordinate 
of  the  bottom  of  the  metal  layer  is  specihed  under  zmin  where  zmin  is  the  distance 
in  microns  from  the  bottom  of  the  substrate  layer  to  the  bottom  of  the  metal  layer. 

You  can  also  dehne  optional  vias  to  make  connections  from  metal  layer  to  metal 
layer.  Note  that  vias  are  totally  optional  and  usually  play  a  minor  role  in  determining 
the  characteristics  of  a  device.  On  the  other  hand,  one  can  create  multi-layer  devices, 
or  symmetric  center-tapped  inductors,  with  many  via  connections  in  series  with  the 
device.  In  these  situations  the  inclusion  of  vias  is  necessary.  Via  layers  are  added  in 
the  ‘via’  sections.  The  top  and  bottom  sections  dehne  the  metal  connectivity  of  the 
via.  Note  that  the  numbers  correspond  to  the  numbers  used  in  the  ‘metal’  sections. 


4.5  Electric 

4.5.1  Introduction 

The  Electric  VLSI  Design  System  is  a  highly  hexible  and  powerful  system  that 
can  handle  many  diherent  types  of  circuit  design  (MOS,  Bipolar,  schematics,  printed 
circuitry,  hardware  description  languages,  etc.)  It  handles  geometry  at  any  angle  (not 
just  Manhattan)  and  can  even  handle  curves. 

Layout  is  done  by  placing  and  wiring  electrical  components.  Although  this  is 
standard  practice  for  schematics,  it  is  unusual  for  chip  layout.  However,  because  of 
this  style  of  design.  Electric  understands  chip  layout  at  a  more  sophisticated  level,  and 
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can  aid  in  design  to  an  unprecedented  degree.  Electric  maintains  circuit  information 
based  on  connectivity.  This  feature  was  made  us  to  incorporate  the  idea  of  local 
reference  terminal  in  Electric. 

Electric  has  many  analysis  tools,  including  design-rule  checking,  simulation,  and 
network  comparison.  Electric  has  many  synthesis  tools,  including  routing,  com¬ 
paction,  silicon  compilation,  PLA  generation,  and  compensation. 

The  user  interface  is  quite  sophisticated  and  runs  on  all  popular  workstations 
(Windows,  Macintosh,  and  UNIX).  It  also  provides  interpretive  languages  (Lisp,  TCL, 
and  Java)  for  advanced  users. 

The  most  interesting  feature  of  the  system  is  its  global  enforcement  of  connectivity 
which  provides  top-down  design  capability  and  ease  of  post-design  modihcations.  The 
features  provided  by  Electric  is  described  in  the  progammers  reference  manual  [7]. 

In  the  next  few  sections  we  will  see  how  a  layout  is  prepared  for  either  EM  simu¬ 
lation  or  in  general  the  process  for  going  from  layout  to  circuit  netlist  simulation. 

4.5.2  Identifying  the  Local  Reference 

A  facet  as  shown  in  Figure  4.10  can  represents  a  local  reference  group  or  a  system 
of  local  reference  groups.  Electric  is  based  on  the  connectivity  approach,  it  stores 
terminal  information  for  the  given  schematic.  This  information  is  used  to  identify  the 
reference  terminals  on  a  facet. 

The  tool  was  modihed  to  automatically  identify  the  local  reference  terminal.  This 
information  is  available  to  the  user  as  a  extension  to  the  GIF  hie  of  the  layout. 

4.5.3  Creating  Exports  in  a  Schematic 

The  user  can  explicitly  name  the  exports  or  terminals  on  a  facet.  This  gives  the 
user  the  advantage  to  identify  local  reference  terminals  in  a  system  of  local  reference 
groups.  This  allows  the  naming  of  nodes  at  points  of  connection  between  schematic 
and  layout  as  shown  in  Figure  4.11. 
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Figure  4.10:  Local  Reference  Terminal. 
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Electric 


Figure  4.11:  Creating  Export. 
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With  the  facility  to  name  ports,  the  user  can  mark  the  internal  and  external  ports 
required  for  EM  simulation  within  the  layout  tool  itself. 


4.5.4  Schematic  Capture 

Up  to  this  stage,  we  have  represented  the  circuit  description  of  the  power  combin¬ 
ing  system  and  identihed  the  local  reference  terminals  on  the  circuit.  Based  on  this 
data,  we  can  generate  a  netlist  describing  the  electrical  circuit  as  shown  in  Figure 
4.12. 
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Figure  4.12:  Creating  spice  netlist. 
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4.6  Integrating  Electric  and  EMPDK 

In  addition  to  the  netlist  we  also  output  two  files,  which  make  the  system  an  EM- 
Aware  system.  These  hies  describe  the  geometry  of  the  distributed  structures  in  the 
layout  and  also  provide  process  parameters  of  the  layers  used  to  layout  the  circuit. 


Mixed 

Layout 

Electric 


Layout 

Parasitic 

Extraction 

_ 

Extractor 

EM-Awcre 
Technology 
file _ 


M/  \l/ 


FREEDA 


Figure  4.13;  Integrating  Electric  and  EMPDK. 


The  integration  of  Electric  and  EMPDK  is  shown  in  Figure  4. 13. The  above  infor¬ 
mation  is  given  as  input  to  fREEDA  to  extract  the  port  voltages  in  the  circuit. 


4.7  fREEDA 

4.7.1  Introduction. 

fREEDA(TM)  is  a  multi-physics  circuit  simulator  under  development  by  a  user 
community  from  universities,  communities  and  laboratories. 

It  uses  state  variables,  local  reference  group  concepts,  automatic  differentiation 
and  hence  captures  multo-physics  and  makes  model  development  straightforward. 
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The  simulation  approach  realizes  the  hrst  new  approach  since  SPICE-like  analyses 
were  developed.  The  approach  allows  the  modeling  of  virtually  any  physical  model 
and  the  generic  model  evaluation  mechanism,  in  which  the  primitive  model  equations 
are  wrapped  in  analysis  specihc  generic  functions,  reduces  the  time  required  not  only 
for  computation  but  also  development. 

fREEDA(TM)  implements  several  types  of  analyses,  all  of  them  state  variable 
based,  including  a  DC  analysis,  time  marching  transient  analysis  (  different  integra¬ 
tion  methods  available),  harmonic  balance  analysis  and  a  unique  wavelet  transient 
analysis. 

It  also  implements  various  linear  devices  from  resistors  to  Foster’s  canonical  port 
representation  and  non-linear  devices  from  diodes  to  three  and  four  terminal  transis¬ 
tors. 

4.7.2  Simulating  structures 

.  Up  to  this  stage,  we  have  complete  description  of  the  layout  under  inspection. 
To  complete  the  suite,  we  required  the  netlist  to  be  in  a  format  that  could  be  read 
into  fREEDA.  Since  Electric  provides  us  with  a  SPICE  netlist,  we  had  to  develop  a 
script  to  automate  the  process  of  conversion.  The  netlist  description  is  given  as  input 
to  fREEDA,  to  generate  the  port  voltages. The  simulation  process  is  shown  in  Figure 
4.14. 


Figure  4.14:  Simulation  of  distributed  structures  under  fREEDA. 
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Chapter  5 


Case  Study 


5.1  Quasi  Grid  Amplifier 

5.1.1  Introduction 

The  earlier  chapters  described  the  approach  towards  an  integrated  design  environ¬ 
ment  for  design  of  spatial  power  combining  systems.  In  this  chapter  various  results 
are  presented.  The  achievement  is  that  the  same  results  are  obtained  with  the  current 
design  environment  compared  to  previous  results  in  [10]. 

In  this  section  various  size  grids  are  analyzed  and  the  intermediate  GIF  data 
is  compared  with  that  in  [10].  The  principle  achievement  is  that  the  GIF  output 
generated  automatically  by  this  work  is  exactly  same  to  that  handgenerated  by  user 
in  [10].  This  indicates  that  the  GIF  output  precisely  captures  the  geometry. 

5.1.2  Unit  Cell 

The  GIF  output  provided  in  [Appendix  A]  is  for  a  unit  cell  shown  in  Fig.  5.1.  Its 
magnitude  and  phase  are  plotted  versus  frequency  in  Fig.  5.2.  The  solid  line  shows 
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Figure  5.1:  Unit  Cell  Dimensions. 

results  obtained  by  simulation  and  the  dashed  line  is  for  the  measurements  taken.  In 
Fig.  5.1  ,  the  geometry  for  a  unit  cell  of  dimension  93.8  mm  x  93.8  mm  with  a  gap 
spacing  of  9.8  mm  and  a  line  width  of  6.35  mm  is  shown. 

5.1.3  5x5  Grid 

The  data  (GIF)  generated  for  this  grid  was  compared  and  found  to  match  with 
the  data  in  [10].  In  [10],  the  design  and  simulation  were  performed  for  the  5x5  grid 
array  placed  in  the  lens  system  at  z  =  0  on  a  dielectric  substrate  with  er  =  2.56  and 
thickness  9.5  mm.  It  may  be  noted  that  the  signihcant  improvement  over  [10]  is  the 
ability  to  automatically  generate  GIF  information  with  local  reference  nodes. 
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FREQUENCY  (GHz) 

Figure  5.2:  Unit  Cell  S21. 

5.2  Microwave  Stub-filter 

5.2.1  Introduction 

Figure  5.2  is  a  microwave  stub  filter.  A  microwave  filter  is  a  two  port  network 
used  to  control  the  frequency  response  at  a  certain  point  in  a  microwave  system  by 
providing  transmission  at  frequencies  within  the  passband  of  the  hlter  and  attenuation 
in  the  stopband  of  the  hlter.  Typical  frequency  response  include  low-pass,  high-pass, 
band-pass  and  band-reject  characteristics.  Applications  can  be  found  in  virtually  any 
type  of  microwave  communication,  radar,  or  test  and  measurement  system. 

There  are  two  methods  to  design  microwave  hlters  as  described  by  Pozar  [12]. 
One  is  the  image  parameter  method,  which  consist  of  a  cascade  of  simpler  two-port 
hlter  sections  to  provide  the  desired  cutoh  frequencies  and  attenuation  characteristics, 
but  do  no  allow  the  specihcation  of  a  frequency  response  over  the  complete  operat¬ 
ing  range.  The  second  is  the  insertion  loss  method,  which  uses  network  synthesis 
techniques  to  design  hlters  with  a  completely  specihed  frequency  response. 
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4.6  mm 


23  mm 


Figure  5.3:  Microwave  stub  filter. 

Both  these  methods  provide  lumped  element  circuits.  For  microwave  applica¬ 
tions,  such  designs  usually  must  be  modihed  to  use  distributed  elements  consisting  of 
transmission  line  sections.  The  Richard’s  transformation  and  the  Kuroda  identities 
provide  this  step. 

Figure  5.3  describes  the  process  or  the  layers  that  makeup  the  stub  hlter. 

5.2.2  Simulation  Results 

In  this  particular  case-study,  we  were  able  to  successfully  identify  distributed 
structures  in  a  layout  as  opposed  to  considering  then  as  simple  wires.  The  individual 
components  in  the  netlist  we  correctly  identihed.  The  user  had  to  manually  enter 
the  analysis  type  required.  Figure  5.4  shows  the  extracted  netlist  of  the  microwave 
stubhlter  under  inspection. 

The  extracted  netlist  were  checked  for  errors  and  simulated  using  fREEDA.  The 
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TOP  COVER 


METAL  LAYER 


BOTTOM  COVER 


Figure  5.4:  Stub  filter  block. 


waveforms  show  the  return  loss  and  the  insertion  loss  of  the  hlter. 
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*transim  netlist 

vsource:vl  1  0  vac  =  1 
.ref  5 

tlinp4:tl  21  0  2  0 

+  zOmag= 75.00  length=5e-3  k=7  tand=.01  fscale=l  .elO 
alpha=l . 

+  nsect  =  20  fopt=35e9 
tlinp4:t2  23  4  5 

+  z0mag=38.53  length=18e-3  k=7  tand=.01  fscale=l  .elO 
alpha=l . 

+  nsect  =  20  fopt=35e9 
tlinp4:t3  30  6  0 

+  zOmag= 75.00  length=10e-3  k=7  tand=.01  fscale=l  .elO 
alpha=l . 

+  nsect  =  20  fopt=35e9 

res:rl  6  32  r=75 
vsource:v2  32  0 
res:r2  1  21  r=75 
res:ropen  4  5  r=10K 


.ac  start  =  4.5e9  stop  =  8e9  njreqs  =  100 

.out  plot  element  "vsource:v2"  0  If  mag  db  In  ''often. out" 
.out  plot  element  "vsource:vl "  0  If  In  "Input.out" 

.end 


Figure  5.5:  Microwave  Stubfilter  extracted  netlist. 
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Frequency  (GHz) 
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.6:  Sll  -  Stub  filter 
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The  above  s-parameters  obtained  from  simulation  were  in  agreement  with  mea¬ 
sured  data  obtained  from  experiments  conducted  in  the  lab. 

5.2.3  Performance  Measures 

Advantages  of  the  presented  tool-flow 

•  Since  all  the  tools  are  GPL  licensed  software,  allows  programmability /customizability 
and  seamless  integration. 

•  Can  be  used  with  additional  libraries  like  MPI (parallel  algorithm)  to  improve 
efficiency 

Accuracy  of  the  simulated  results 

It  has  been  shown  that  the  automatically  generated  GIF  hie  matches  with  the  hand 
generated,  published  by  Patwardhan[10].  With  the  incorporation  of  the  local  reference 
groups,  we  make  sure  that  effects  like  time  delay,  substrate  coupling  are  simulated 
accurately.  The  current  distribution  values  and  the  s-parameters  extracted  from  the 
full-wave  simulator  for  simple  structures  were  in  agreement  with  the  measured  data. 


5.3  Summary 

The  above  mentioned  case  studies  tested  the  validity  of  the  individual  tools  in 
the  framework.  The  GIF  output  generated  in  the  hrst  case-study  was  in  agreement 
with  the  hand  generated  GIF  hie.  The  simulated  results  for  the  structure  were  in 
agreement  with  the  measured  data.  In  the  second  case-study,  the  tool  was  tested  for 
its  accuracy  in  capturing  distributed  structures.  The  tool  was  able  to  identify  most 
of  the  structures  that  were  laid  out. 
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Chapter  6 


Conclusion 


An  integrated  electromagnetic  and  circuit  simulation  environment  capable  of  sim¬ 
ulating  high  speed  circuits  has  been  implemented.  The  analysis  incorporates  surface 
modes  and  full  nonlinear  effects.  A  4  cell  grid  ampliher  was  used  to  validate  the 
simulation  tools  with  measurements.  The  package  consists  of  a  method  of  moments 
full-wave  electro-magnetic  simulator,  along  with  a  general  nonlinear  circuit  simula¬ 
tion  program.  The  programs  and  the  interfaces  between  have  been  dehned  and  doc¬ 
umented.  The  simulation  environment  allows  grid  ampliher  structures  to  examined 
from  a  wholistic  point  of  view.  The  effect  of  edge  effects,  mutual  coupling,  and  device 
variations  can  now  be  considered.  Given  a  generic  layout  hie,  we  can  attach  it  with 
any  process.  The  simulation  environment  will  allow  what-if  studies  to  be  conducted, 
limiting  the  amount  of  empirical  design  required  Inherent  problems  in  the  grid  design 
such  as  edge  elements  can  now  be  considered.  The  ehects  of  device  variations  can  also 
be  examined.  Design  ideas  can  now  be  implemented  on  the  computer  hrst,  increasing 
the  overall  chance  of  success 
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Chapter  7 


Future  Research 


There  are  many  areas  that  could  be  pursued  further  with  regard  to  this  simula¬ 
tion  environment.  The  programs  in  the  hybrid  co-simulation  environment  need  to 
be  combined  together  in  a  more  coherent  fashion.  There  is  a  lot  of  scope  to  imple¬ 
ment  to  EM  analysis  using  FDTD  techniques.  Support  for  alternate  incident  held 
distributions  should  also  be  added.  Currently  only  uniformly  distributed  helds  are 
considered.  Other  held  distributions,  such  as  gaussian  distributions  should  also  be 
implemented.  Oscillator  and  stability  analysis  also  need  to  be  included  in  the  tool 
set.  This  would  require  the  inclusion  of  an  oscillator  analysis  in  the  circuit  simulator. 
The  oscillation  analysis  should  also  include  a  means  to  test  the  stability  of  a  grid 
ampliher  system.  Another  improvement  in  the  circuit  simulator  is  the  ability  to  use 
subcircuits.  Subcircuits  are  currently  supported  for  passive  elements.  This  needs  to 
be  extended  to  include  nonlinear  elements  also.  The  use  of  subcircuits  becomes  very 
important  as  the  size  of  the  structures  examined  increases.  The  chances  of  errors  in 
the  netlists  is  greatly  reduced  and  the  speed  of  making  changes  is  greatly  increased. 
Along  with  improvements  in  the  circuit  simulator,  the  active  device  models  also  need 
to  be  improved.  A  strategy  for  developing  nonlinear  behavior  models  is  essential.  The 
results  from  the  circuit  simulator  will  in  most  cases  be  limited  by  the  accuracy  of  the 
active  device  model.  This  is  especially  true  under  large  signal  operation.  Stability 
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analysis  will  require  an  active  device  model  which  is  accurate  over  a  large  frequency 
range.  The  same  requirements  will  hold  for  designing  high  efficiency  power  ampliher 
systems  where  harmonic  loading  is  important. 
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Chapter  8 


Appendix 


8.1  Appendix  A 

8.1.1  GIF  file  format 

Caltech  Intermediate  Format  (GIF)  is  a  recent  form  for  the  description  of  inte¬ 
grated  circuits.  Created  by  the  university  community,  GIF  has  provided  a  common 
database  structure  for  the  integration  of  many  research  tools.  GIF  provides  a  limited 
set  of  graphics  primitives  that  are  useful  for  describing  the  two-dimensional  shapes 
on  the  different  layers  of  a  chip.  The  format  allows  hierarchical  description,  which 
makes  the  representation  concise.  In  addition,  it  is  a  terse  but  human-readable  text 
format.  GIF  is  therefore  a  concise  and  powerful  descriptive  form  for  VLSI  geometry. 

Each  statement  in  GIF  consists  of  a  keyword  or  letter  followed  by  parameters  and 
terminated  with  a  semicolon.  Spaces  must  separate  the  parameters  but  there  are  no 
restrictions  on  the  number  of  statements  per  line  or  of  the  particular  columns  of  any 
field.  Comments  can  be  inserted  anywhere  by  enclosing  them  in  parenthesis. 

There  are  only  a  few  GIF  statements  and  they  fall  into  one  of  two  categories: 
geometry  or  control.  The  geometry  statements  are:  LAYER  to  switch  mask  layers. 
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BOX  to  draw  a  rectangle,  WIRE  to  draw  a  path,  ROUNDFLASH  to  draw  a  circle, 
POLYGON  to  draw  an  arbitrary  figure,  and  CALL  to  draw  a  subroutine  of  other 
geometry  statements.  The  control  statements  are  DS  to  start  the  definition  of  a 
subroutine,  DF  to  finish  the  definition  of  a  subroutine,  DD  to  delete  the  definition  of 
subroutines,  0  through  9  to  include  additional  user-specified  information,  and  END 
to  terminate  a  GIF  file.  All  of  these  keywords  are  usually  abbreviated  to  one  or  two 
letters  that  are  unique. 

Label  Command 

94  Input  10,20  Type’Tnput”  Reference=”  Ground” ; 

This  command  is  another  common  extension  and  is  used  to  mark  points  with  a 
name.  The  parameters  are  the  name  of  the  label,  and  the  x  and  y  coordinates  of  the 
point. 

Extensions  to  GIF  can  be  done  with  the  numeric  statements  0  through  9.  Al¬ 
though  not  officially  part  of  GIF,  certain  conventions  have  evolved  for  the  use  of 
these  extensions 

Sample  GIF  File 

This  is  a  sample  GIF  file  modeling  a  dipole  antenna.  The  units  are  in  millimeters. 
The  name  of  the  layout  is  ‘dipole’.  It  consists  of  two  boxes  on  a  metall  type  layer. 

The  centers  of  two  cells  are  located  at  the  coordinates  (-100;  0)  and  (100;  0), 
and  both  cell  sizes  are  200  mm  x  200  mm  (square  boxes).  A  circuit  port  (See  chapter 
3  for  port  description)  is  specified  at  the  coordinate  (0;  0). 

DS  1  1  1; 

9  Dipole; 

L  metall; 

B  200  200  -100,0; 

B  200  200  100,0; 

94  Ground  0,0  Type=” Ground” ; 

DF; 

C  1; 

E 
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8.2  Appendix  B 

8.2.1  EMPDK  User  Guide 

8.2.2  Introduction 

EMPDK  prepares  a  layout  file  for  EM  simulation.  Some  of  the  basic  requirements 
are  the  material  properties  of  the  layers  in  the  layout.  The  output  of  EMPDK  is  in  a 
human-readable  form.  This  hie  can  be  imported  to  other  layout  tools  like  Cadence, 
LinkCad  etc. 

8.2.3  Creating  a  new  technology  file 

The  user  must  have  information  regarding  the  process  and  the  layers  that  go  into 
the  layout.  To  create  a  new  hie,  the  user  has  to  select  the  option  ‘New’  from  the  ‘File’ 
Menu. 

After  providing  the  ‘hlename’  the  user  will  be  prompted  for  information  regarding 
the  layers  in  the  process.  This  information  has  be  entered  starting  with  the  substrate 
layer  working  your  way  upwards. 

Once  the  user  is  done  entering  the  information  for  each  of  the  layers  in  the  process, 
he/she  can  select  the  ‘Done’  button  in  the  dialog  box.  This  will  automatically  generate 
the  hie  for  the  user. 

More  information  regarding  the  layer  types  supported  and  how  to  specify  dimen¬ 
sions  are  given  in  Appendix  C. 

8.2.4  Editing  an  existing  technology  file 

To  edit  an  existing  hie,  the  user  must  select  the  option  ‘Open’  from  the  ‘File’ 
Menu.  Once  the  user  provides  the  tool  with  the  layout  hlename,  the  tool  displays  a 
dialog  box  with  the  layers  used  for  the  layout. 
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The  user  must  have  process  information  regarding  the  layers  generated  by  the 
layout  tool.  This  information  can  be  entered  in  the  dialog  box.  Once  the  user  is 
done,  the  process  can  be  completed  by  pressing  the  ‘Done’  button  in  the  dialog  box. 
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8.3  Appendix  C 


8.3.1  Technology  File  Structure 


//TECHNOLOGY  FILE 

II  METAL  LAYERS  :  2 
II DIELEGTRIO  LAYERS:  1 


technology  cmos 
UNITS  microns 
metal 


//name 

zmln 

thickness  reslstlvlty(ohms,cm)  desc 

layer  0 

0,00 

0,35  0,2  //G  substrate 

layer  1 
endmetal 

0,7  0,2 

1 ,56E-4  II M 

dielectric 

//name 

zmln 

thickness  dielectric  constant 

layer  2 

0,35 

0,35  3,9 

enddlelectrlc 

endtechnology 


Figure  8.1:  Sample  technology  file 


The  technology  hie  allows  the  user  to  specify  the  process  related  parameters.  The 
characteristic  information  obtained  from  the  technology  is  used 

•  to  specify  the  thickness  of  the  layers 

•  obtain  the  ordering  of  the  layers 

•  to  obtain  the  electrical  parameters  like  resistivity,  dielectric  constant,  etc. 

The  structure  of  the  hie  is  very  simple  and  human  readable.  The  hrst  line  is  a 
comment  line  which  gives  us  the  version  nnmber  of  the  technology  hie.  The  next  two 
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lines  tell  the  simulator  the  number  of  metal  and  dielectric  layers  in  this  particular 
process.  The  description  of  the  process  starts  with  the  process  name  specihed  by  the 
’technology’  tag.  Then  we  provide  the  simulator  with  the  units  in  which  the  related 
distance  information  is  provided.  By  default  it  is  microns.  Then  we  specify  the  metal 
layers  between  the  tags  ’metal’  and  ’endmetal’.  The  information  is  specihed  in  the 
bottom-up  fashion,  i.e,  we  provide  the  substrate  and  then  the  layer  closest  to  it  and 
work  our  way  upwards.  In  the  similar  fashion  we  provide  information  of  other  layers 
like  ’dielectric’  and  ’via’.  The  end  of  the  hie  is  denoted  by  the  ’endtechnology’  marker. 
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